An ab-initio many-body method for electronic structure calculations of solids. 
II. Unscreened Hartree-Fock treatment for the 3d systems Fe, Co, Ni and Cu 
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The ab-initio many-body method suggested in the preceding paper is applied to the 3d tran- 
sition metals Fe, Co, Ni, and Cu. We use a linearized muffin-tin orbital calculation to determine 
Bloch functions for the Hartree one-particle Hamiltonian, and from these obtain maximally localized 
Wannier functions. Within this Wannier basis all relevant one-particle and two-particle Coulomb 
matrix elements are calculated. The resulting second-quantized many-body Hamiltonian with ab- 
initio parameters is studied within the simplest many-body approximation, namely the unscreened, 
selfconsistent, Hartree-Fock approximation (HFA). We present these HFA results, which we believe 
are the first to have been done for crystalline 3d transition metals, and compare them with those 
obtained from the standard local (spin) density approximation (LSDA) within density functional 
theory (DFT). Although the d-bands sit considerably lower within HFA than within L(S)DA, the 
exchange splitting and magnetic moments for ferromagnetic Fe, Co, and Ni are only slightly larger 
in HFA than what is obtained experimentally or within LSDA. The HFA total energies are lower 
than the corresponding L(S)DA calculations. 

PACS numbers: 71.10.Fd, 71.15.AP, 71.15.Mb, 71.20.Bo ,71.45.Gm, 75.10.Lp 
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I. INTRODUCTION 

In the preceding paperi we suggested a new method for 
ab-initio electronic-structure calculations of solids. The 
main steps of this procedure are: 

1. Perform a conventional, self-consistent, band- 
structure calculation for an effective one-particle 
Hamiltonian, namely, the Hartree Hamiltonian, to 
obtain a suitable basis set of Bloch functions. 

2. By taking into account only a finite number J of 
bands one chooses a truncated one-particle Hilbert 
space. The Marzari-Vanderbilt^ algorithm is then 
used to construct a maximally localized set of Wan- 
nier functions, which span the same truncated one- 
particle Hilbert space. 

3. All one-particle (tight-binding) and two-particle 
(Coulomb) matrix elements of the Hamiltonian 
within this Wannier function basis are calculated. 

4. The resulting electronic many-body Hamiltonian in 
second quantization with parameters determined 
from first principles is studied within standard 
many-body approximations for lattice electron sys- 
tems. 

In this paper we apply this scheme to the 3d transition 
metals Fe, Co, Ni, and Cu. We use the "linear mufhn- 
tin orbital" (LMTO) method within the "atomic-sphere 
approximation" (ASA)2i to perform the band-structure 
calculation for the Hartree Hamiltonian in first quanti- 
zation. The direct Coulomb matrix elements of the max- 
imally localized Wannier basis are rather large, about 
20 eV in magnitude. We then use the simplest possible 



many body approximation, the Hartree-Fock approxima- 
tion (HFA) , to study the second-quantized multiple-band 
Hamiltonian. Our results are compared with those ob- 
tained from a standard LSDA calculation^"*'- Although 
the 3d-bands and the 4s-band overlap in the L(S)DA ap- 
proximation, our unscreened HFA calculations give 3d- 
bands that lie considerably lower (between 10 and 20 eV) 
than the 4s-band. The HFA correctly predicts ferromag- 
netism for the ferromagnetic metals Fe, Co, and Ni and 
no magnetism for Cu, but with a much larger exchange 
splitting between majority and minority 3d bands than 
obtained within LSDA and with a slightly larger mag- 
netic moment per site than obtained experimentally or 
within LSDA. On the other hand, the total energy is 
lower in HFA than in LSDA. The LSDA results for met- 
als are probably more reliable than our new HFA results, 
which lack important screening and correlation effects. In 
order for our method to go beyond LSDA we would need 
to use better many body methods than the (unscreened) 
HFA, which should be possible within our scheme. 

Our purpose in presenting the HFA results is to demon- 
strate that the new ab-initio many-body method that we 
have proposed is feasible and can be applied to practical 
calculations of materials. In addition, because HFA is the 
assumed standard starting point before adding complex 
many-body correlations, and because it is the simplest 
many-body approximation (or the best mean-field one- 
particle approximation), it is useful to know what the 
HFA predicts for the 3d transition metals, which have 
been so heavily studied by other techniques. Compar- 
isons of HFA results with experiments, other (higher or- 
der) theories, or established standard methods such as 
the L(S)DA should demonstrate the effects of correlation 
on electronic properties in d-electron systems. 
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To the best of our knowledge we do not know of any 
published HFA results (band structure, density of states, 
magnetism, magnetic moment, total energy, etc.) for 
the 3d ferromagnets Fe, Co and Ni, unless it was im- 
plicitly applied to these materials for schemes like the 
local ansatz^:, where HFA results serve as an input to 
higher order calculations . The is not surprising since 
the HFA has, from very early on, been viewed as a poor 
approximation for metals. For example, when applied 
to the homogeneous electron gas (as the simplest model 
of an infinite metallic system), the HFA has well-known 
Fermi edge singularitieai^. These lead, in particular, to 
a vanishing density of states (DOS) at the Fermi en- 
ergy, which is, of course, unphysical. This unphysical 
feature usually prevails in actual HFA-calculations for 
real metals^, though sometimes this singularity is hard 
to see in actual HFA-results'^°. In our calculations the 
non-locality is handled through the calculation of expec- 
tation values (matrix elements of the density matrix), 
which makes HF calculations as easy as Hartree calcula- 
tions. Furthermore, because of our localized Wannier ba- 
sis, we only keep on-site and a few inter-site Coulomb and 
exchange matrix elements. Hence our calculations have 
an effective short-ranged Coulomb interaction. Although 
longer-range Coulomb matrix elements are small in our 
calculations, which is why we truncate them, it is possi- 
ble that if all of them were kept to infinite distances that 
they could add up to give Fermi edge singularities (which 
are due to the long-ranged nature of the bare Coulomb 
interaction) and other standard anomalies. Correlation 
or screening would quickly kill these effects. 

The paper is organized as follows. In Section^we de- 
scribe the LMTO-Hartree calculations and the Hartree- 
results for the band structure and density of states, re- 
spectively. Section IIIII describes some results obtained 
for the maximally localized Wannier functions obtained 
within the Marzari-Vanderbilt algorithm^, in particular 
their localization properties. Results for the matrix el- 
ements, in particular the direct Coulomb and exchange 
matrix elements are given in Section IIVI we also com- 
pare these results with calculations of the Slater inte- 
grals. The application of the (unscreened) HFA to the 
multiband many-body Hamiltonian in second quantiza- 
tion is the subject of Section For an interpretation 
of the results we compare the numerical HFA results ob- 
tained for the crystal with previous atomic HFA results 
and with numerical and analytical results for a simplified 
local (atomic or zero band width) model in Section IVII 
A comparison with the more standard LSDA-results fol- 
lows in Section rVlII before the paper closes with a short 
discussion. 



II. LMTO HARTREE CALCULATION 

For the four materials of interest (Fe, Co, Ni, Cu) we 
performed a selfconsistent Hartree band-structure calcu- 
lation. Besides the nuclear charge we used the (experi- 
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FIG. 1: Total density of states (of both spin directions) 
around the Fermi level (chosen as zero of the energy axis 
Ef = 0) obtained from the selfconsistent Hartree calculation 
for Fe. 

mentally) known results for the lattice structure (bcc for 
Fe, fee otherwise; Co should actually be hexagonal) and 
for the lattice constant as input. For the band-structure 
calculation we used the LMTO method^-H within the 
atomic sphere approximation (ASA). The radius of the 
(overlapping) muffin-tin spheres (the Wigner-Seitz radius 
S) is determined by the condition that the sphere volume 
equals the volume of the unit cell. This yields the fol- 
lowing values for the Wigner-Seitz radius: S = 2.662ao 
for Fe, S = 2.621ao for Co,^ S" = 2.602ao for Ni and 
S = 2.669ao for Cu (Ref. i). Within the muffin-tin 
spheres the potential and wave functions are expanded 
in spherical harmonics with a cutoff 1^3.^ — 2, i.e., s, p, 
and d-orbitals are included. 

Within the LMTO-ASA the eigenfunctions, i.e. the 
Bloch wave ftmctions (inside a muffin-tin sphere), are 
given in terms of the solution to the radial Schrodinger 
equation (f),yi (r) to some fixed energies -B^; and its energy 
derivative ipviir)'- 

*nk(r)-^(0.,(^M2'' + '^./Wi?f)i^L(f) , (1) 
L 

where the Yl = denote the (complex) spherical har- 
monics and n is the band index. 

Figs, n - 0] show the nonmagnetic (spin-degenerate) 
Hartree results for the density of states (DOS). One ob- 
serves that the five narrow 3d-bands are similar in energy 
and therefore hybridize with the 4sp-bands. For Fe, Co, 
and Ni the Fermi level Ep falls in the midst of the par- 
tially filled 3d-bands, and hence there is a large DOS at 
the Fermi level. For Cu Ep is above the filled 3d bands, 
where the states have more 4sp-like character and the 
DOS is small (free electron like). 

The Bloch functions obtained from this Hartree calcu- 
lation are used as the basis of a truncated Hilbert space. 
In the band calculation, there is already a natural trun- 
cation due to the spherical harmonics expansion cutoff 
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FIG. 2: The same as in Fig. [Hfor Co. 
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FIG. 4: The same as in Fig. [T]for Cu. 
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FIG. 3: The same as in Fig. [T]for Ni. 



atom per unit cell. The connection between the Bloch 
functions in the different muffin-tin spheres is given by 
the standard Bloch relation 



■0„k(r + R) 



ikR 



^nk(r) 



(2) 



Hence the Bloch function in a single muffin-tin sphere 
determines the function for the whole crystal. However, 
this simple relation holds for Bloch functions only. Any 
other function (to be represented as linear combination of 
Bloch functions) can, nevertheless, be decomposed into 
its contributions within the different muffin-tins accord- 
ing to 



$„(r) =^$„(R,;r-R,) 



(3) 



where $q(R; r) = for |r| > S. 



III. MAXIMALLY LOCALIZED WANNIER 
FUNCTIONS 

We now turn to the construction and description of the 
Wannier functions and will present results for their main 
properties for the four materials under investigation. Be- 
cause 



WR«(r) = u)o„(r - R) , 



(4) 



it is sufficient to investigate the Wannier functions only 
for one arbitrarily chosen lattice vector 0: w„(r) = 
woji(r). This Wannier function has contributions not 
only in the muffin-tin sphere around but also in other 
muffin-tin spheres, and we can use the decomposition 
(PI into the contributions from the different muffin-tin 
spheres 



Wn{r) =^w„(R;r-R) 



(5) 



R 



After application of the Marzari-Vanderbilt algorithm, 
which is briefly summarized in the preceding papesi and 
in more detail in Ref. 0, the new set of bands that are 
used to calculate the maximally localized Wannier func- 
tions can no longer be classified by pure angular mo- 
mentum quantum numbers. The Wannier functions are 
rather admixtures having different L-contributions (3d, 
4s, 4p etc.). But, since the original Bloch functions from 
which the Wannier functions are constructed were given 
in terms of a spherical harmonics expansion, the new 
Wannier functions (and their contribution in each indi- 
vidual muffin-tin sphere) can also be decomposed into 
these spherical harmonics contributions 



(R;r) = ^{0.z(r)A^-f0.,(r)i3f"}ri(f) . (6) 

L 



One can then calculate the weight of the contributions to 
the Wannier function (centered at 0) within the different 
muffin-tin spheres 



(wn|wn)R= / dh\wnir)\'^ ^ / d^r | w„ (R; r) | ^ , 

JR Jo 



(7) 
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n 





1 


2 


3 


4 


5 


6 


7 




.9761 


.9765 


.9596 


.9800 


.9773 


.8754 


.8731 


.8763 




.0019 


.0018 


.0081 


.0019 


.0017 


.2224 


.2381 


.2265 




.0955 


.0726 


.1797 


.0611 


.0728 


.5480 


.5509 


.5347 




.9026 


.9256 


.8121 


.9370 


.9255 


.2295 


.2110 


.2388 



TABLE I: Some properties of the lowest eight maximally 
localized Wannier functions of Fe. 



and one can also decompose this into the different l- 
contributions according to: 



{WnW, 



i)r 



E 



E {l^^f + (8) 



For the 3d-system iron these quantities are tabulated in 
Table m The first line is the weight {wn\wn)o in the cen- 
ter muffin-tin. Between 88 and 98% of the total weight of 
the Wannier functions is to be found already within the 
center muffin-tin; this shows how well localized our Wan- 
nier functions are with the lowest five functions having 
values of more than 95%. Rows 2-4 in this table indicate 
the different Z-contribution or Z-character of the Wannier 
functions. One sees that the optimally localized Wan- 
nier functions arc not pure within their Z-character, but 
the lowest five Wannier functions (0-4) still have mainly 
I = 2 (3d) character. Higher band-index states (which 
are slightly less well localized according to row 1) are ad- 
mixtures that have mainly I = 1 (4p) character (about 
50 %), but also a considerable amount of / = (4s) and 
I = 2 (3d) character. Corresponding results for the other 
3d-systems Co, Ni, and Cu are similar and, therefore, not 
repeated here. 





22.42 
20.90 
20.10 
20.96 
20.86 
14.16 
13.32 
13.96 
13.50 



20.90 
23.04 
19.95 
21.55 
21.53 
14.07 
13.54 
13.58 
14.15 



20.10 
19.95 
20.77 
20.05 
19.83 
12.95 
13.46 
13.37 
13.22 



20.96 
21.55 
20.05 
23.27 
21.67 
13.46 
14.05 
13.98 
13.98 



20.86 
21.53 
19.83 
21.67 
22.99 
13.71 
13.28 
14.25 
14.12 



14.16 
14.07 
12.95 
13.46 
13.71 
13.67 
9.45 
9.58 
9.64 



6 

13.32 
13.54 
13.46 
14.05 
13.28 
9.45 
13.52 
9.27 
9.50 



7 

13.96 
13.58 
13.37 
13.98 
14.25 
9.58 
9.27 
13.75 
9.65 



13.50 
14.15 
13.22 
13.98 
14.12 
9.64 
9.50 
9.65 
13.81 



Jnm 





1 


2 


3 


4 


5 


6 


7 


8 





22.42 


0.84 


0.61 


0.75 


0.99 


0.86 


0.73 


0.81 


0.42 


1 


0.84 


23.04 


0.77 


0.88 


0.84 


0.70 


0.51 


0.48 


0.86 


2 


0.61 


0.77 


20.77 


0.88 


0.70 


0.96 


0.93 


0.92 


0.60 


3 


0.75 


0.88 


0.88 


23.27 


0.82 


0.33 


0.78 


0.64 


0.69 


4 


0.99 


0.84 


0.70 


0.82 


22.99 


0.52 


0.46 


0.75 


0.83 


5 


0.86 


0.70 


0.96 


0.33 


0.52 


13.67 


0.58 


0.56 


0.57 


6 


0.73 


0.51 


0.93 


0.78 


0.46 


0.58 


13.52 


0.45 


0.56 


7 


0.81 


0.48 


0.92 


0.64 


0.75 


0.56 


0.45 


13.75 


0.55 


8 


0.42 


0.86 


0.60 


0.69 


0.83 


0.57 


0.56 


0.55 


13.81 


TABLE IL 


On-site direct and exchange Coulomb matrix el- 


ements between Wannier functions for Fe 


AU 


energ 


les are 



eV's. 



Wannier functions are explicitly given as linear combi- 
nations of products of spherical harmonics and a radial 
wave function. The expansion 
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- J2 m^') YKir) (11) 
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IV. ONE PARTICLE AND COULOMB MATRIX 
ELEMENTS 

From the optimally localized Wannier functions we cal- 
culate the one-particle matrix elements 

ti2 = J d?Y wI{y) (^-^V2 + V[V)^ W2{V) (9) 

and the Coulomb matrix elements of the Hamiltonian 

Wi2M = j (i'vd'v' wl{v)w*{v') w:i{v')wi{v) . 

(10) 

Here we use the abbreviated notation 1 to mean Ririi 
and 2 to mean for R2n2, etc. We used two different 
numerical algorithms to calculate these Coulomb matrix 
elements, namely the FFT-algorithm briefly described in 
the previous paper and a spherical expansion algorithm 
described in some detail in Ref. 12. The latter method 
makes use of the fact that (in each muffin-tin sphere) the 



[K = {fc, m}) makes it possible to express the on- 
site Coulomb integrals as one-dimensional integrals over 
products of the radial functions and the (tabulated) 
Gaunt coefficients. The results obtained by this al- 
gorithm and by the independent FFT-algorithm agree 
within the numerical errors (of at most 1%). 

Results for the on-site direct and exchange Coulomb 
matrix elements between the optimally localized Wan- 
nier functions are given in Table ^ for iron (Fe) . The 
direct Coulomb integrals Unm = Wnm.mn between the 
Wannier states with the lowest five band indices (n, m G 
{0, . . . ,4}), which according to the discussion and table 
in the preceding section have mainly 3d-character, are 
rather large, up to 23 eV for Fe. Within the 3d-like 
bands the interband direct Coulomb matrix elements are 
of the same magnitude as the intraband matrix elements. 
The matrix elements between 3d-states and 4sp-states 
are considerably smaller, of the magnitude of 13 - 14 eV. 
For electrons in 4sp-states (n, m G {5, . . . , 9}) the direct 
intraband Coulomb matrix elements are again of the or- 
der of 13 - 14 eV, but the interband matrix elements are 
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U 


J 


tNN 


tNNN 


Fe 


21.1 


.81 


.59 


.24 


Co 


22.6 


.87 


.55 


.10 


Ni 


22.6 


.88 


.75 


.11 


Cu 


24.5 


.94 


.80 


.12 



TABLE III: Averaged on-site Coulomb, exchange, nearest 
neighbor and next nearest neighbor hopping matrix elements 
for the 4 3d-systems; energies are in eV. 





pO p2 pA 


Fe (crystal) 
Fe (atom [14]) 


21.62 9.61 5.91 
23.76 10.96 6.81 


Co (crystal) 
Co (atom [Ml) 


23.18 10.31 6.34 
25.15 11.58 7.20 


Ni (crystal) 
Ni (atom [14]) 


24.69 11.00 6.77 
26.53 12.20 7.58 


Cu (crystal) 
Cu (atom [14]) 


26.27 11.72 7.23 
27.90 12.82 7.96 



TABLE IV: Slater integrals (in eV) for the 3d-systems 
Fe, Co, Ni, Cu as obtained by our calculations and within an 
earlier atomic calculatioij—. 



slightly smaller, about 9 eV. The exchange matrix ele- 
ments Jnm = Wnm.nrn are always much smaller, usually 
less than 1 eV (for n ^ m). Again the corresponding 
results for the other 3d-systems investigated (Co, Ni and 
Cu) are very similar. 

For the 5 states with predominant 3d-character we 
have calculated the averages of the on-site direct and 
exchange Coulomb matrix elements 



U 



mm 

_ 1 

'J = "^Q / ^ ^^min' mm' i 



(12) 
(13) 



as well as the averages of the absolute values of the near- 
est neighbor (NN) and next nearest neighbor (NNN) hop- 
ping matrix elements 



NN{N) 



1 ^ 



(14) 



The results obtained thereby for the 4 transition metals 
under consideration are shown in Table lTTll The U- values 
vary between 21 eV for Fe and 25 eV for Cu, the J-values 
are smaller than 1 eV and the hopping matrix elements 
are of the magnitude 0.5 - 0.7 eV for nearest-neighbor 
(NN) and 0.1 - 0.2 eV for next-nearest-neighbor (NNN), 
and further on decrease with increasing distance. 
We have also evaluated the Slater integralsi^: 



dr / dr' r'^ |i?/=2(r)|^ 



where Ri=2{f) is a radial (atomic) d-wave function (ob- 
tained by solving the Schrodinger equation for a radial 
symmetric potential, for instance). Note that only the 
three integrals , and are required to determine 
all the Coulomb d-matrix elements. Using the radial d- 
wave function obtained from the Hartree calculation we 
obtain the following values for the Slater integrals of the 
four 3d-systems: F° = 21.62 eV for Fe, 23.18 eV for Co, 
24.69 eV for Ni, and 26.27 eV for Cu. This means, the 
Slater integrals F^ are rather good estimates of our (av- 
eraged) Coulomb matrix elements. These values are also 
in agreement with older results obtained in calculations 
for 3d-atomsi4. In Table HVl we show our i^*^ -values for 
the four 3d-crystals and compare them with correspond- 
ing atomic calculations from Ref. ^3- Obviously, there 
is fairly good agreement between these atomic and our 
results. 



V. UNSCREENED HARTREE-FOCK 
APPROXIMATION 

After we have determined the matrix elements within 
our restricted basis set of 9 maximally localized Wannier 
functions (per site and spin), we have a Hamiltonian in 
second quantization of the form 

H = ^ti2c\^C2a + ^ W^12,34cL4^/C3„'C4^ (16) 
12cr 1234crcr' 

for which all the matrix elements are known from first 
principles. The simplest approximation one can now ap- 
ply is the HFA, which replaces the many-body Hamilto- 
nian by the effective one-particle Hamiltonian 



12(T 

^HF _ ^Hart , ^Fock 



(17) 



with S-!". = Sf/- + Ef2°f (18) 



34(t' 



(15) 



Here the expectation values {c[^C2a) have to be deter- 
mined selfconsistently for the HF Hamiltonian (|17() . Note 
that the Fock (exchange) term is spin (a) dependent and 
may, therefore, give rise to magnetic solutions. 

The Hartree-Fock results for the four materials of in- 
terest are shown in Figs. [511121 We show the effective 
HF band structure and its density of states (DOS). In 
our HFA calculations there are no singularities (or a van- 
ishing DOS) at the Fermi level since we start from a lo- 
calized description and consider the Coulomb matrix ele- 
ments only locally (on-site and for a few neighbor shells). 
Therefore, we implicitely truncate the Coulomb interac- 
tion in real space and in practice work with an effective 
short-ranged interaction. Within HFA the main part of 
the 3d-bands lies between 18 and 22 eV below the Fermi 
level and is separated from the 4sp-bands. We find mag- 
netism in HFA for Fe, Co, and Ni in agreement with 
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FIG. 5: (Unscreened) Hartree-Fock band-structure of Fe; 
the full line shows the majority (spin up) , the dashed line the 
minority spin component. 



FIG. 7: The same as Fig. 1^1 for Go. 
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Energy relative to Fermi energy [eV] 

FIG. 6: Density of states (per spin direction) for Fe within 
the (unscreened) Hartree-Fock approximation; the full line 
shows the majority (spin up), the dashed line the minority 
spin contribution. 



experiment. The five majority spin d-bands are about 
20 eV below the Fermi energy and are completely filled. 
But the partially filled minority d-bands have two (for 
Fe), three (for Co), and four (for Ni) filled bands be- 
tween -18 and -15 eV, and the rest are around and above 
the Fermi level. This results in magnetic moments per 
atom (in units of the Bohr magneton /i^/atom) of 2.9 for 
iron, 1.9 for cobalt, and 0.76 for nickel. For copper no 
magnetism and exchange splitting of the 3d-bands is ob- 
tained, but the (spin degenerate) 3d-bands are at about 
22 eV below the Fermi level and separated from the 4sp- 
bands. If we compare these results with the results of the 
simple Hartree approximation shown in Figs, ^^we see 
that the exchange term has two effects: it produces an 
exchange splitting and the possibility of magnetic solu- 
tions, and it draws the 3d-bands energetically down by an 
amount of about 20 eV. Compared with experiment the 
HFA overestimates magnetism and leads to overly large 
values for the magnetic moment per site (about 10-30% 
too large; experiment gives magnetic moments of about 
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FIG. 8: The same as Fig. El for Go. 



2.2 for iron, 1.7 for Co, and 0.6 for Ni). This is consis- 
tent with Heisenberg or Ising model studies where the 
mean-field approximation HFA also has the tendency to 
overestimate magnetism and magnetic solutions. 

However, the reason why the 3d-bands lie so far be- 
low the Fermi level and the 4sp-band in HFA has noth- 
ing to do with the existence and overestimation of mag- 
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FIG. 9: The same as Fig. Elfor Ni. 
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FIG. 11: The same as Fig. Elfor Cu. 



netism. This can be seen already from the non-magnetic 
system Cu, for which the (fully occupied) 3d-bands also 
lie at about 22 eV below the Fermi level (see Figs. ^2 
and I12|) . To demonstrate this also for a system with 
a partially filled 3d-band we have done a non-magnetic 
Hartree-Fock-calculation for Co (by forcing equal occu- 
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FIG. 12: The same as Fig. IHlfor Cu. 



FIG. 14: Total density of states (of both degenerate spin 
directions) for Co obtained within the non-magnetic HFA- 
solution. 

pation for both spin directions). The results for the band 
structure and the DOS are shown in Figs. I13I14I We ob- 
serve again that the main part of the 3d-bands arc well 
below the 4s-bands and Fermi level; note the hybridiza- 
tion gap caused by the unoccupied 3d-bands above the 
Fermi level. 



VI. COMPARISON WITH ATOMIC 
HARTREE-FOCK RESULTS 

We have seen in the previous section that one effect of 
the HFA calculation, when compared with the Hartree 
calculation, is the shift of the 3d-bands down (about 
20 eV below the Fermi level and about 8-10 eV below 
the bottom of the 4sp-band). This shift of the d-bands 
is about the same energy as the Coulomb matrix ele- 
ments U, and roughly agrees earlier atomic Hartree-Fock 
calculationsi^ii^, where the 3d-states are also about 10 
eV below the 4s-states. 

Because the inter-site hopping matrix elements in Ta- 
ble are much smaller than the U- values one may con- 
sider an expansion in t/U, with the zeroth order ap- 
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proximation to completely neglect hopping. Doing this, 
we have performed a quasi-atomic HFA calculation for 
Co, by including only the on-site one-particle and two- 
particle (Coulomb) matrix elements. The results are 
summarized in Fig. ^1 the degeneracy of the different 
levels is also indicated. In the paramagnetic case, we find 
that the 3d-bands are below the 4s-bands (at the Fermi 
level) by about 6 to 7 eV, which is in rough agreement 
with the earlier atomic HFA results^'*-^^. The splitting 
between the occupied and unoccupied 3d-states is about 
23 eV, which is the on-site U for Co. Magnetic HFA 
solutions are also found in the atomic limit for Co, as 
shown in the right panel of Fig. 1151 The majority-spin 
3d states {T2g and Eg) are now completely filled and en- 
ergetically lie lower than the corresponding non-magnetic 
HFA-states. But only the (3-fold degenerate) T2g-states 
of the minority-spin electrons are filled whereas the Eg- 
states of the minority electrons are empty (and now even 
26 eV above the occupied d-states) . The additional ener- 
getical shifts between the occupied 3d-states in the para- 
magnetic and ferromagnetic atomic HFA solution are due 
to the exchange matrix elements J. 

This behavior can qualitatively be understood within 
the framework of the following simple, analytically solv- 
able model. Similar to the numerical HFA-results pre- 
sented and discused above, we neglect all intersite one- 
particle (hopping) and interaction matrix elements. Fur- 
thermore, we assume that we have diagonalized the one- 
particle Hamiltonian, taking into account only the atomic 
3d-levels and assuming that the on-site one-particle diag- 
onal matrix elements e, the Coulomb matrix elements U, 
and the exchange matrix elements J are equal, i.e., that 
the 3d-levels are degenerate in the atomic limit with no 
crystal-field effects. Then the atomic part of the many- 
body Hamiltonian can be written as 



+ 2 H 4o-c]<T'Cm'Cj<T (19) 

where i,j G {0, ...,4} denote the 5 (degenerate) 3d- 
states. 

The standard Hartrcc-Fock decoupling leads to 
H = J2{e + C/[E(4'C..')-(cU.)] 

ia jcr' 

- •/E^sV^J-))clc.^ • (20) 

Here we have assumed that the Hartree-Fock Hamil- 
tonian has the same symmetry as the uncorrelated 
Hamiltonian, and hence off-diagonal expectation values 
(c]^/Cicr) for («cr) ^ (jcr') vanish. From this equation it is 
clear that the HF Hamiltonian can be written in terms 
of an effective one-particle energy 

^ = E^?/M^- (21) 

where 

jo' j^i 

(22) 

In the simple Hartree approximation (HA) the exchange 
decouplings are neglected, which means that all the de- 
coupling terms with the negative sign would not occur. 
Therefore, the corresponding Hartree one-particle ener- 
gies are given by 

ef/ = £ + C/E(sVc,.') • (23) 

JO-' 

Comparing this result with the Hartree-Fock one-particle 
energies, we find that the HF occupied levels are shifted 
downwards by an amount of 

relative to the Hartree levels. Momentarily setting J = 0, 
we see that for N occupied levels the Hartree approxima- 
tion gives the one-particle energies 

ef/ = e + NU (25) 

whereas the HFA yields 

efJ^ = e+{N~l)U . (26) 

The occupied Hartree-Fock one-particle energies are 
lower than the corresponding Hartree one-particle ener- 
gies by [/, which is a consequence of the artificial and 
unphysical self-interaction still present in the Hartree ap- 
proximation that is exactly cancelled in Hartree-Fock. 
This also explains why the Hartree-Fock bands are 
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shifted downwards from the Hartree bands by an en- 
ergy of the amount U. One also sees from this simple 
atomic-limit Hartree-Fock model that the energy differ- 
ence between the highest occupied and the lowest un- 
occupied effective Hartree-Fock one-particle energies is 
again essentially U, which is once more in agreement 
with our numerical results for the crystal and for the 
atom (cf. Fig ll5|l . Note that we have ignored Usd inter- 
actions, which cause an additional shift of d-bands below 
the s-bands by about an additional 10 eV in the full HFA 
calculations. 

Taking into account the exchange interaction J again 
and denoting by Ng- the number of occupied states with 
spin (7 (i.e. N = + Ni) one obtains in HFA 



.HFA 



= s+iN -1)U -iNg-l)J 



Then the total energy in HFA is given by 
NiN ^Y^N^N^- 



Etnt - Ne 



-J 



(27) 



(28) 



For the total energy we have added the necessary cor- 
rection term to the sum of the occupied energy levels 
(much like the double counting term that shows up in 
band-structure calculations). Now for partially filled 3d- 
shells the occupation of the different spin directions may 
be different. Denoting M = N-^ — Ni we obtain for the 
total energy 

N(N ~ 1) N'^ 
Etot =Ne + '-U - 



^^^^^ ^ m 



The magnetic {M ^ 0) total energy is lower than the 
nonmagnetic (consistent with Hund's rules). 

Take once more Co with 8 3d-clectrons. The param- 
agnetic (nonmagnetic) state has the occupations N \^ — 
A'^i = 4 {N=8 and M=0). For this configuration (corre- 
sponding to the left panel in Fig. I15|l one obtains 



8e + 28U-12J 



(30) 



The Hund's rule magnetic solution has 3d-statcs of one 
spin-direction completely filled, i.e., = 5 and = 3 
{N=8 and M=2). This gives 



E, 



(M) 



^ 8e-f 28C/- 13J 



(31) 



Therefore, the magnetic configuation (with a magnetic 
moment of 2 for the atom) is energetically more favorable 
by J . Note also the exchange splitting in the occupied 
energy eigenvalues 



(32) 



and that our model would predict the unoccupied minor- 
ity spin E2g state to be J7 -I- J higher in energy than the 
corresponding occupied majority spin state. 

The simple model in this section differs from the results 
shown in Fig. 1151 in that we have replaced the full matrix 
of U and J by scalar values for d-states only (ignoring s-d 
interactions, for example). However, it captures all of the 
important physics without attempting to be completely 
quantitative. 
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FIG. 17: The same as Fig. EHlfor Co. 

COMPARISON WITH LSDA RESULTS 



For comparison with the HFA results described in Sec- 
tion we have also performed a standard local spin- 
density approximation (LSDA), ab-initio, band-structure 
calculation with the LMTO-ASA method. The results for 
the DOS are shown in Figs. 1 161 1 191 For the magnetic sys- 
tems Fe, Co, and Ni one obtains also an exchange split- 
ting and the prediction of magnetic solutions with mag- 
netic moments of 2.18 for iron, 1.58 for Co, and 0.58 for 
Ni, which are in better agreement with experiment than 
the HFA results. This is consistent with the expectation 
that L(S)DA is tuned to accurately predict ground-state 
properties. 

The energy spectra of the bands (DOS) are quite dif- 
ferent from the HFA. For example, the 3d-bands now 
fall into the same energy region as the 4sp-bands, i.e., 
the LSDA-results are not so different from the Hartree- 
results. This means that the exchange-correlation energy 
leads only to a small shift of the 3d-bands downwards by 
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FIG. 19: Total density of states (of both degenerate 
spin directions) for Cu within the local density approxima- 
tion(LSDA) 



at most a few eV and a smaller exchange splitting (also 
of the magnitude of 1 eV) in the case of magnetic solu- 
tions. On the other hand, in the LSD A calculations the 
self- interaction terms are not completely canceled, i.e., an 
(unrealistic) self-interaction is included, which may lead 
to 3d bands that lie energetically too high, as discussed 
for the atomic limit in the previous section. 

To sec the effect of correlations within LSDA, we have 
also performed an exchange-only calculation for Co. The 
result is shown in Fig. [501 Obviously the (majority) d- 
bands lie lower than the ones in the full LSDA-calculation 
shown in Fig. El But the shift is of the magnitude of 1 
eV only, i.e., very minimal when compared with the large 
drop in the full HFA. On the other hand, this exchange- 
only LSDA-result also does not contain self-interaction 
corrections, which are responsible for the large shift 
downwards of the d-bands. Nevertheless, the LSDA re- 
sult indicates that a possible effect of correlations is to 
shift the 3d-bands up relative to exchange-only calcu- 
lations, and hence one would expect a similar effect if 
correlations could be added to the full HFA-calculations. 
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3d transition metals Fe through Cu. 



One can also calculate the total energy in the Hartree, 
HFA, and LSDA approximations. The results obtained 
for the four materials of interest are shown in Fig. |^ 
We see that the total energy is always significantly lower 
in HFA than in the Hartree approximation, which is ex- 
pected because the HFA minimizes the total energy. The 
HFA total energy is also lower than the L(S)DA, and 
the LSDA result is lower than the simple Hartree re- 
sult. Because of the unknown approximations that go 
into constructing L(S)DA, it is hard to guess ahead of 
time that this would be the case. However, it is well 
known that the L(S)DA approximation produces a bad 
total exchange-correlation energy; the reason why such 
good agreement with experiment is found is that relative 
exchange-correlation energies are nonetheless reasonably 
accurately calculated. 
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VIII. DISCUSSION AND CONCLUSION 

In this paper we have presented tlie results of (un- 
screened) HFA calculations for the 3d transition met- 
als Fe, Co, Ni and Cu. We obtain magnetic solutions 
for Fe, Co, and Ni with (slightly) too large magnetic 
moments when compared to experimental or LSDA re- 
sults. The occupied HFA 3d-bands lie about 20 eV be- 
low the Fermi level (and the Hartree result), which is 
also the magnitude of the splitting between occupied and 
unoccupied 3d-bands and of the magnitude of the on- 
site Coulomb matrix element (the "Hubbard" U). This 
downwards shift of the HFA 3d-bands compared to the 
Hartree- and L(S)DA-3d-bands can be understood as due 
to the self-interaction correction of HFA. 

One may argue that these results are not surprising 
and an artifact of using the unscreened HFA. Our ab- 
initio calculation of the direct Coulomb matrix elements 
yields large values of the magnitude of 20 eV. There- 
fore the Coulomb energies are large for these materials 
and HFA can be considered to be an approximation for 
the selfenergy which is correct only in linear order in 
the Coulomb interaction. But for these large values of 
the U-terms HFA is certainly not sufhcient but one has 
to apply better many-body approximations. One should 
apply systematic extensions of HFA, which within the 
standard perturbational approach can be represented by 
(a resummation of an infinite series of) Feynman dia- 
grams, or one can try to apply the recently so successful 
non-perturbational many-body schemes like "dynamical 
mean field theory" (DMFT)i& or variational (Gutzwiller) 
approaches''^^. The simplest standard diagram series are 
the bubble diagrams leading essentially to the "random 
phase approximation" (RPA). This means just a renor- 
malization of the interaction line, i.e. the pure "naked" 
Coulomb interaction has to be replaced by a "dressed" 
interaction. Or in other words, the exchange (Fock) con- 



tribution has not to be calculated with the bare Coulomb 
matrix elements but with screened Coulomb matrix el- 
ements. Probably the non-perturbational schemes like 
DMFT are also only applicable for screened Coulomb 
matrix elements. 

To summarize, the new many-body ab-initio approach 
proposed in Ref. Ill is applicable to real materials. Maxi- 
mally localized Wannier functions and their one- and two- 
particle matrix elements can be calculated from first prin- 
ciples so as to obtain a many-body multi-band Hamil- 
tonian in second quantization. Standard methods of 
many-body theory can then be applied to this Hamil- 
tonian. The simplest approximation (HFA) is unreli- 
able, because the large Coulomb matrix elements are un- 
screened. Therefore, in the future better many-body ap- 
proximations should be used. Within the standard Feyn- 
man diagram approach the most straightforward next 
step consists in a summation of bubble diagrams leading 
to a renormalized (screened) Coulomb interaction. This 
requires calculating the exchange contribution not with 
the bare but with a screened Coulomb interaction. To 
take into account the effects of screening one has to cal- 
culate the charge susceptibility and the (static) dielectric 
constant, which can be done within a generalized Lind- 
hard theory, for instance. 
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